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ABSTRACT 

SN HFF14Tom is a Type la Supernova (SN) discovered at z = 1.3457 ± 0.0001 behind the galaxy cluster Abell 
2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 ± 0.15 magnitudes brighter 
than unlensed Type la SNe at similar redshift, implying a lensing magnification of /r 0 bs = 2.03 ± 0.29. This observed 
magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 
lens models, 13 of which were generated before the SN magnification was known, qualifying as pure “blind tests”. 
The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with 
the measured fj, to within ler. However, there is a subtle systematic bias: the significant disagreements all involve 
models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or 
methodological explanation to account for it. We do find that model accuracy can be improved to some extent with 
stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. 
In addition to testing model accuracies as we have done here, Type la SN magnifications could also be used as inputs 
for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional 
strong- and weak-lensing information is unavailable. 

Subject headings: supernovae: general, supernovae: individual: HFF14Tom, galaxies: clusters: gen¬ 
eral, galaxies: clusters: individual: Abell 2744, gravitational lensing: strong, grav¬ 
itational lensing: weak 
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1. INTRODUCTION 

Galaxy clusters can be used as cosmic telescopes to 
magnify distant background objects through gravita¬ 
tional lensing, which can substantially increase the reach 
of deep imaging surveys. The lensing magnification en¬ 
ables the study of objects that would otherwise be unob- 
servable because they are either intrinsic ally faint (e.g. 


distant (e.g. Franx et al. 1997 

Ellis er al. 2001 

1 Hu 

et al.12002; Kneib et al.j2004; Ric 

rard et al. 2006, 

2008 


Bouwens et al.|2009t|Maizy et al" 2010 Z heng et al. 2012. 

Coe et al.|2013| 'Bouwens et al.||20T4 Zitrin et al.||2014 ). 

Background galaxies are also syatiaily magnified , aliow- 


early universe with resolutions of ^100 pc (e.g. Stark 

et al. 20081 

Jones et al.| 2010; Yuan et al. 2011; Wuyts 

eTaT 2014 

jivermore et al. 2015|. 


dow onto the transient sky through an appropriately 
cadenced imaging survey. The flux magnification from 
strong-lensing clusters is especially valua ble for the st udy 


of z > 1.5 supernovae (SNe) 

(e-g- 

Kovner & Paczynski 

1988; Kolatt & Bartelmann 1998; Sullivan et al. 2000f 

Saini et al. 2000; Gunnarsson & Coobar 2003; Coobar 

et al. 2009) Postman et al. 

2012), which are still ex- 


Riess et al. 2001, 2007; 

Suzuki et al. 2012 

Rodney et al. 

2012 Rubin et al.j[2013 

Jones et al. 2013 

)• 


imental setup: instead of using strong-lensing clusters to 
study distant SN e, we ca n use the SNe as tools for exam¬ 
ining the lenses (Riehm ~et al.| [2011). The most valuable 


transients for testing and improving cluster lens models 
would be stro ngly tensed, SNe that are resolved in to mul¬ 
tiple images (Holz 2001 Oguri Sz Kawano||2003 1. Tran¬ 


sients that are lensed into multiple images can also be¬ 
come cosmological tools, as the measurement of time de¬ 
lays between the images can provide cosm ograp hic infor¬ 


mation to constrain the Hubble para meter (Refsdal 19641 
and other cosmological parameters (Linder ' 201 If . We 


have recentl y observed the fir st example of a multiply- 


imaged SN ( Kelly et al.||2015 ). We expect to detect the 


reappearance of this object (called “SN Refsdal”) within 


the next year 1 

Oguri 2015 

Sharon & Johnson 2015| 

Diego 

et al. 

2015), d 

elivering a precise test ol lens mod 

el pre- 


rently very unlikely (Li et al. 2012), the y will become 


much more common m the nex t decade (Coe & Mous- 
takas 2009 Dobke et al. 2009), and may be developed 


into an imp ortant new cosmological tool (Oguri 2010 


Linder 


2011 ). 


In addition to time delays from multiply-imaged SNe, 
we can also put cluster mass models to the test with the 


magnified but not multiply-imag 

3d. 

Patel et al. 

(2014 

hereafter P14) and Nordin et al. 

(201 

_4j) presented inde- 


2 are securely classified as Type la SNe all found in 
the Cluster Lensing and Supernova survey with Hubble 


(CLASH, P hPostman, HST Program ID 12068, Postman 
et al.][2012 ). Both groups demonstrated that these stan- 


dard candles can be used to provide accurate and precise 
measurements of the true absolute magnification along a 


random sight line through the cluster. Although in these 
cases the SNe were used to test the cluster mass models, 
one could in principle incorporate the measured magni¬ 
fications of Type la SNe into the cluster as additional 
model constraints. In that role, Type la SNe have the 
particular value that they can be found anywhere in the 
cluster field. Thus, they can deliver model constraints in 
regions of “middle distance” from the cluster core, where 
both strong- and weak-lensing constraints are unavail¬ 
able. Moreover, given enough time, multiple background 
Type la SNe could be measured behind the same cluster, 
each providing a new magnification constraint. 

One of the key values in observing standard candles 
behind gravitational lenses is i n addressing the problem 
of the mass-sheet degeneracy (Falco et al. 1985 Schnei¬ 
der & Seitz 1995). This degeneracy arises because one 
can introduce into a lens model an unassociated sheet 
of uniform mass in front of or behind the lens, without 
disturbing the primary observable quantities. For exam¬ 
ple, take a lens model with a given surface mass den¬ 
sity k, and then transform the surface mass density to 
k' = (1— X)k+X for any arbitrary value A. Both the k and 
k! models will produce exactly the same values for all po- 
sitio nal and shear constraints from strong and weak lens¬ 
ing ( Seitz fe Schneider||1997 ). When lensed background 
sources are available across a wide range of redshifts (as 
is the case for the Abell 2744 cluster discussed here), it 
should in principle be possible to break this degeneracy 
(Seitz & Schneider |1997 Bradac et al.||2004 ). However, 
there are mor e com plex versions of positional constraint 
degeneracies (|Liesenborgs fe de Rijckc 2012; |Schneider| 
& Sluse 2014). Such degeneracies do not extend to the 
absolute magnification of a background source’s flux and 
size. In the case of the simple mass sheet degeneracy de¬ 
scribed a bove, the ma gnific ation scales as p oc (1 — A)~ 2 
(see e.g. Bartelmann 2010). Therefore, one can break 
these fundamental degeneracies with an absolute mea¬ 
surem ent of magnification from a standard candle (Holz 
20011 or a standard ruler (Sonnenfeld et al.||2011). 

In Section [2] we present the discovery and follow-up 
observations of SN HFF14Tom at z = 1.3457, discovered 
behind the galaxy cluster Abell 2744. Section [3] exam¬ 
ines the SN host galaxy. Sections 0] and [5] describe the 
spectroscopy and photometry of tnis SN7 leading to a 
classification of the object as a normal Type la SN. In 
Section [6] we make a direct measurement of the magni¬ 
fication of this source due to gravitational lensing. Sec¬ 
tion [7] discusses the tension between our magnification 
measurement and the lens models. Our conclusions are 
summarized in Section [8j along with a discussion of fu¬ 
ture prospects. 

2. DISCOVERY, FOLLOW-UP, AND DATA PROCESSING 

SN HFF14Tom was discovered in Hubble Space Tele¬ 
scope ( HST) observations with the Advanced Camera 
for Surveys (ACS) in the F606W and F814W bands (V 
and I), collected on UT 2014 May 15 as part of the 
Hubble Frontier Fields (HFF) survey (PLJ.Lotz, HST- 
PID:13495) The HFF program is a 3-year Direc¬ 
tor’s discretionary initiative that is collecting 140 or¬ 
bits of HST imaging (roughly 340 ksec) on six mas¬ 
sive galaxy clusters, plus 6 accompanying parallel fields. 

36 http:/ /www. stsci.edu/hst/campaigns/frontier-fields 
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Figure 1. SN HFF14Tom in the Abell 2744 field. The left panel shows a UV/Optical/IR color composite image constructed from all 
available HST imaging of the Abell 2744 cluster field. The inset panels on the right show F814W imaging of the immediate vicinity of SN 
HFF14Tom, approximately 40 // from the center of the cluster. The top panel shows the template image, combining all data prior to the 
SN appearance. Labeled ellipses mark the nearest galaxies and their spectroscopic redshift constraints, with the most likely host galaxy 
marked in blue, and a background galaxy in yellow. The bottom panel is constructed from all HFF F814W imaging taken while the SN was 
detectable, and marks the SN location with an arrow. (Left panel image credit: NASA, ESA, and J. Lotz, M. Mountain, A. Koekemoer, 
and the HFF Team (STScI)) 


Table 1 

J2000 Coordinates of HFF14Tom, host, and cluster. 


Object 

R.A. 

(h:m:s) 

Decl. 

(d:m:s) 

R.A. 

(deg) 

Decl. 

(deg) 

HFF14Tom 

00:14:17.87 

-30:23:59.7 

3.574458 

-30.399917 

Host galaxy 

00:14:17.88 

-30:24:00.6 

3.574483 

-30.400175 

Abell 2744 a 

00:14:21.20 

-30:23:50.1 

3.588333 

-30.397250 


a Coordinates of the HFF field center, approximately at the center 
of the cluster. 


Each field is observed in 3 optical bands (ACS F435W, 
F606W and F814W) and 4 infrared (IR) bands (WFC3- 
IR F105W, F125W, F140W, and F160W), although the 
optical and IR imaging campaigns are separated by 
months. Abell 2744 was the first cluster observed, with 
IR imaging spanning 2013 October-November, and opti¬ 
cal imaging from 2014 May-July. A composite image of 
the HFF data showing the SN is presented in Figure [l] 
and the locations of the cluster center, the SN, and the 
presumed host galaxy are given in Table [l] The SN de¬ 
tection was made in difference images constructed using 
template imaging of Abell 2744 from HST+ACS obser¬ 
vations taken in 2009 (PLDupke, HST-PID:11689). 

Upon discovery, HST target-of-opportunity obser¬ 
vations were triggered from the FrontierSN program 
(PPRodney, HST-PID:13386), which aims to discover 
and follow transient sources in the HFF cluster and paral¬ 
lel fields. The FrontierSN observations provided WFC3- 
IR imaging as well as spectroscopy of the SN itself using 
the ACS G800L grism, supplementing the rapid-cadence 
optical imaging from HST+ACS already being provided 


by the HFF program. The last detections in the IR 
F105W and F140W bands came from the direct-imaging 
component of the GLASS program. Difference images 
for the IR follow-up data were generated using templates 
constructed from the HFF WFC3-IR imaging campaign, 
which concluded in November, 2013. 

All of the imaging data were processed using the 
sndrizpipe pipeline^ 7 ] a custom data reduction pack¬ 
age in Python that employs the DrizzlePac tools Jrom 
the S pace T elescope Science Institute (STScI) (Fruchter 
et al.[2010). Photometry was collected using the PyPhot 
software package a pure-Python implementation of 
the p hotometry algo rithms from the IDL AstroLib pack¬ 
age (Landsman 1993), which in tur n are based on the 
DAOPHU'l program (Stetson 1987). For the IR bands 
we used point spread function (PSF) fitting on the differ¬ 
ence images, and in the ACS optical bands we collected 
photometry with a ty/3 aperture. Table [2lpresents the list 
of observations, along with measured photometry from 
all available imaging data. 


3. HOST GALAXY 

The most probable host galaxy for SN HFF14Tom is 
a faint and diffuse galaxy immediately to the south-east 
of the SN location. With photometry of the host galaxy 
collected from the template images, we fit the spectral en¬ 
ergy distribution (SED) using t he BPZ code - a Bayesian 
photometric redshift estimator (Bemtez|2000). From the 


https://github.com/srodney/sndrizpipe vl.2 

DOI:10.528l/zenodo.lU731 

38 https://github.com/djonesl040/PyPhot 
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Table 2 

HFF14Tom Observations and Photometry 


Obs. Date 
(MJD) 

Camera 

Filter 
or grism 

Exp. Time 
(sec) 

Flux 

(counts/sec) 

Flux Err 
(counts/sec) 

AB Mag a 

Mag Err 

AB Zero Point 

AZP b 

(Vega-AB) 

56820.06 

ACS 

F435W 

5083 

-0.027 

0.053 

27.66 


25.665 

-0.102 

56821.85 

ACS 

F435W 

5083 

0.105 

0.053 

28.11 

0.55 

25.665 

-0.102 

56823.77 

ACS 

F435W 

5083 

0.022 

0.053 

29.80 

2.59 

25.665 

-0.102 

56824.97 

ACS 

F435W 

5083 

0.021 

0.053 

29.85 

2.72 

25.665 

-0.102 

56828.68 

ACS 

F435W 

5083 

-0.148 

0.053 

27.65 


25.665 

-0.102 

56830.87 

ACS 

F435W 

5083 

0.100 

0.054 

28.16 

0.58 

25.665 

-0.102 

56832.86 

ACS 

F435W 

5083 

-0.080 

0.053 

27.66 


25.665 

-0.102 

56833.86 

ACS 

F435W 

5083 

-0.002 

0.053 

27.67 


25.665 

-0.102 

56839.50 

ACS 

F435W 

5083 

-0.022 

0.052 

27.68 


25.665 

-0.102 

56792.06 

ACS 

F606W 

5046 

0.363 

0.083 

27.59 

0.25 

26.493 

-0.086 

56792.98 

ACS 

F606W 

3586 

0.692 

0.095 

26.89 

0.15 

26.493 

-0.086 

56797.10 

ACS 

F606W 

4977 

0.968 

0.087 

26.53 

0.10 

26.493 

-0.086 

56800.08 

ACS 

F606W 

4977 

0.844 

0.085 

26.68 

0.11 

26.493 

-0.086 

56804.99 

ACS 

F606W 

5046 

0.977 

0.086 

26.52 

0.10 

26.493 

-0.086 

56792.99 

ACS 

F814W 

3652 

1.639 

0.104 

25.41 

0.07 

25.947 

-0.424 

56797.11 

ACS 

F814W 

4904 

3.376 

0.141 

24.63 

0.05 

25.947 

-0.424 

56798.95 

ACS 

F814W 

5046 

3.951 

0.156 

24.46 

0.04 

25.947 

-0.424 

56800.10 

ACS 

F814W 

4904 

3.854 

0.155 

24.48 

0.04 

25.947 

-0.424 

56801.89 

ACS 

F814W 

10092 

4.102 

0.153 

24.41 

0.04 

25.947 

-0.424 

56802.95 

ACS 

F814W 

10092 

4.325 

0.160 

24.36 

0.04 

25.947 

-0.424 

56803.93 

ACS 

F814W 

15138 

4.402 

0.160 

24.34 

0.04 

25.947 

-0.424 

56804.08 

ACS 

F814W 

5046 

4.658 

0.178 

24.28 

0.04 

25.947 

-0.424 

56812.08 

ACS 

F814W 

637 

4.705 

0.258 

24.27 

0.06 

25.947 

-0.424 

56815.93 

ACS 

F814W 

446 

4.026 

0.285 

24.43 

0.08 

25.947 

-0.424 

56820.07 

ACS 

F814W 

5044 

3.508 

0.142 

24.58 

0.04 

25.947 

-0.424 

56821.87 

ACS 

F814W 

5044 

3.541 

0.144 

24.57 

0.04 

25.947 

-0.424 

56823.79 

ACS 

F814W 

5044 

2.876 

0.124 

24.80 

0.05 

25.947 

-0.424 

56824.99 

ACS 

F814W 

5044 

3.060 

0.129 

24.73 

0.05 

25.947 

-0.424 

56828.70 

ACS 

F814W 

5044 

2.777 

0.121 

24.84 

0.05 

25.947 

-0.424 

56830.89 

ACS 

F814W 

5044 

2.395 

0.111 

25.00 

0.05 

25.947 

-0.424 

56832.88 

ACS 

F814W 

5044 

2.331 

0.108 

25.03 

0.05 

25.947 

-0.424 

56833.88 

ACS 

F814W 

5044 

2.389 

0.111 

25.00 

0.05 

25.947 

-0.424 

56839.52 

ACS 

F814W 

5044 

1.673 

0.093 

25.39 

0.06 

25.947 

-0.424 

56833.14 

WFC3-IR 

F105W 

756 

7.504 

0.239 

24.08 

0.03 

26.269 

-0.645 

56841.82 

WFC3-IR 

F105W 

756 

5.822 

0.208 

24.36 

0.04 

26.269 

-0.645 

56850.06 

WFC3-IR 

F105W 

756 

3.952 

0.207 

24.78 

0.06 

26.269 

-0.645 

56860.62 

WFC3-IR 

F105W 

1159 

2.899 

0.167 

25.11 

0.06 

26.269 

-0.645 

56886.63 

WFC3-IR 

F105W 

1159 

1.216 

0.147 

26.06 

0.13 

26.269 

-0.645 

56891.67 

WFC3-IR 

F105W 

356 

0.971 

0.324 

26.30 

0.36 

26.269 

-0.645 

56893.20 

WFC3-IR 

F105W 

712 

0.954 

0.242 

26.32 

0.28 

26.269 

-0.645 

56954.64 

WFC3-IR 

F105W 

356 

0.521 

0.388 

26.98 

0.81 

26.269 

-0.645 

56817.08 

WFC3-IR 

F125W 

1206 

8.459 

0.191 

23.91 

0.02 

26.230 

-0.901 

56833.15 

WFC3-IR 

F125W 

756 

7.753 

0.255 

24.01 

0.04 

26.230 

-0.901 

56841.83 

WFC3-IR 

F125W 

806 

6.015 

0.227 

24.28 

0.04 

26.230 

-0.901 

56850.07 

WFC3-IR 

F125W 

806 

4.343 

0.224 

24.64 

0.06 

26.230 

-0.901 

56891.86 

WFC3-IR 

F140W 

712 

2.578 

0.344 

25.42 

0.14 

26.452 

-1.076 

56893.06 

WFC3-IR 

F140W 

712 

3.026 

0.363 

25.25 

0.13 

26.452 

-1.076 

56955.58 

WFC3-IR 

F140W 

1424 

1.218 

0.269 

26.24 

0.24 

26.452 

-1.076 

56817.09 

WFC3-IR 

F160W 

1206 

4.831 

0.263 

24.24 

0.06 

25.946 

-1.251 

56833.21 

WFC3-IR 

F160W 

756 

3.965 

0.241 

24.45 

0.07 

25.946 

-1.251 

56841.84 

WFC3-IR 

F160W 

756 

3.011 

0.234 

24.75 

0.08 

25.946 

-1.251 

56850.08 

WFC3-IR 

F160W 

756 

2.744 

0.223 

24.85 

0.09 

25.946 

-1.251 

56860.67 

WFC3-IR 

F160W 

1159 

1.895 

0.177 

25.25 

0.10 

25.946 

-1.251 

56886.64 

WFC3-IR 

F160W 

1159 

1.677 

0.191 

25.38 

0.12 

25.946 

-1.251 

56812.0 

ACS 

G800L 

3490 







56815.7 

ACS 

G800L 

6086 








a For non-positive flux values we report the magnitude as a 3 -a upper limit 
b Zero point difference: the magnitude shift for conversion from AB to Vega magnitude units. 
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BPZ analysis, we found the host to be most likely an ac¬ 
tively star-forming galaxy at a redshift of z = 1.5 ± 0.2. 
This photo-z was subsequently refined to a spectroscopic 
redshift of z = 1.3457 ± 0.0001, based on optical spec¬ 
troscopy of the host (Mahler et al. in prep.) that shows 
two significant (> 10cr) emission lines at 8740.2 and 
8746.7 Angstroms, consistent with the [Oil] AA 3726- 
3729 A doublet. 

The next nearest galaxy detected in HST imaging is 
2.2"northeast of the SN position. It has a redshift of 
z = 1.742 determined from a spectrum taken with the 
G141 grism of the HST WFC3-IR camera, collected as 
part of the Grism Lens-Amplifi ed S urvey fro m Space 
(GLASS, P LTreu, PID:13459 |Treu et al.pOlj |Schmidt 
et al. 2014] ). As we will see in Sections [|and|5) b oth the 
spectroscopic and photometric data from the SN itself 
are consistent with the redshift of the fainter galaxy at 
0 = 1.3457, and incompatible with ^ = 1.742 from this 
brighter galaxy. This means that the latter galaxy is a 
background object and therefore has no impact on the 
SN magnification. 

The galaxy identified as the host is not close enough to 
the cluster core to necessarily be multiply-imaged, but it 
is still possible that the host galaxy is one of the outer 
images of a multiple image system. In such a case, and 
given the position of the SN host, one would expect that 
another image of the galaxy would be present at a similar 
brightness, and would therefore be detectable in HST 
imaging. To date, no plausible candidate for a counter¬ 
image has been identified. 

To measure the stellar mass of the host galaxy we use 


Eq. 8 of Taylor et al. (2011), which relates the rest-frame 
(g-i) color and i-bam 


uminosity to the total stellar mass. 
To derive these values, we fixed the redshift at z = 1.3457 
and repeated the SED fitting using BPZ. From the best- 
fit SED we extracted rest-frame optical magnitudes, and 
corrected them for lensing using a magnification factor of 
/z = 2.0 - a value that we will derive from the SN itself in 
Section [g] From this we determine the host galaxy mass 
to be 10™ M 0 . 

4. SPECTROSCOPY 

A spectrum of SN HFF14Tom was collected with the 
ACS G800L grism on 2014 June 4 and 7, when the SN 
was within 3 observer-frame days of the observed peak 
brightness in the F814W band. The observations listed 
at the bottom of Table [2] - used 5 HST orbits from the 
FrontierSN program for a total spectroscopic exposure 
time of ~10 ksec. The grism data were processed and the 
target spectrum was extracted using a custom pipeline 
I Brammer et al. 2012), which was developed for the 3D- 
HST program (PI: Van Dokkum; PID: 12177, 12328) and 
also used by the GLASS team. 

Figure [2] shows the composite 1-D ACS grism spec¬ 
trum, combining all available G800L exposures, overlaid 
with SN model fits that will be described below. The 
spectrum is largely free of contamination, because the 
orientation was chosen to avoid nearby bright sources 
and the host galaxy is diffuse and optically faint. Thus, 
the SN spectral features can be unambiguously identi¬ 
fied, most notably the red slope of the continuum and a 
prominent absorption feature at ~8700A. 

Spectra of Type I SNe (including all sub-classes la, 
lb and Ic) are dominated by broad absorption features, 
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Figure 2. Redshift and age determination from spectral template 
matching to the the SN HFF14Tom maximum light spectrum. The 
y axis plots flux in arbitrary units, and the x axis marks wavelength 
in A with the observer-frame on the bottom and rest-frame on the 
top. The HFF14Tom spectrum observed with the HST ACS G800L 
grism is shown in black, overlaid with model fits derived from a 
library of Type la templates that have extended rest-frame UV 
coverage. The top two panels show the best matching templates 
when using a smooth 3rd-order polynomial to warp the shape of the 
template pseudo-continuum, with the redshift fixed at z=1.3457 
and then allowed to float as a free parameter. The lower three 
panels show matches found when the templates are not warped, 
both with and without fixing the redshift. The bottom panel shows 
the best match in this set, although at z = 0.98 it is inconsistent 
with the host galaxy redshift prior and the light curve. 
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which cannot in general be use d to d irectl y ext ract a 
spectroscopic redshift (see e.g. Filippenko 1997). As 
described below, we fit template spectra to the SN 
HFF14Tom data in two steps. First we determine a spec¬ 
tral classification and get a preliminary estimate of the 
redshift and age - using the SuperNova IDcntification 
(SNID) software (Blondin & Tonry||2007). Second, we 
refine the redshift and age measurement using a custom 
Type la spectral template matching program. 

4.1. Classification with SNID 

The SNID program is designed to estimate the type, 
redshift, and age of a SN spectrum through cross¬ 
correlation matching wi th a library of templat e spectra, 
using the algorithm of Tonry & Davis (1979). To ac¬ 
count for possible distortions in the broad shape of the 
SN pseudo-continuum due to dust or instrumental cal¬ 
ibration effects, SNID divides each SED by a smooth 
cubic spline fit. This effectively removes the shape of the 
SN pseudo-continuum to leave behind a flat SED super¬ 
imposed with spectral absorption and emission features. 
It is these features which drive the cross-correlation fit, 
so the SNID approach is insensitive to the overall color 
of the SED. We used v2.0 of the SNID template library, 
which includes template SEDs covering all Type la and 
core collapse sub-classes, and has recently been updated 
with co rrections a nd improvements to the Type Ib/c tem¬ 
plates ( Liu fe Modjaz|[2014 ). 

In Si\ ID the goodness of fit is evaluated primarily 
through the rlap parameter, which measures the degree 
of wavelength overlap and the strength of the cross¬ 
correlation peak. Typically, an rlap value > 5 is required 
to be considered an acceptable match. 

To match the SN HFFMTom spectrum we use con¬ 
servative constraints on age and redshift: limiting the 
age to ±5 rest-frame days from peak brightness and 
0.8 < z < 1.8, consistent with the SN light curve and 
the two plausible host galaxies. With these constraints 
we find that the only acceptable match is a normal Type 
la SN near 2 = 1.3. The best match has rlap= 8.7, using 
the normal Type la SN 2005cf at z = 1.35 and age=-2.2 
rest-frame days before peak. In contrast, the best non-la 
matches all have rlap< 2.5. 

Using SNID we can find an acceptable core collapse SN 
match only when we remove all age and redshift con¬ 
straints. In this case the best non-la match is the Type 
Ic SN 1997ef, which delivers rlap= 6.8 at z = 0.51 and 
age=47.3 rest-frame days past peak (71 observer-frame 
days). This is not as good a fit as the best Type la mod¬ 
els, is at odds with the host galaxy redshift prior, and 
is strongly disfavored bvthe shape and colors of the SN 
light curve (see Section [5]). 

From the preceding analysis, we conclude that 
HFFMTom is a Type la SN at z ~ 1.35. At this redshift, 
the absorption at ~8700A corresponds to the blended 
Ca II H&K and Si II A3858 features. Generally referred 
to as the Ca H&K feature, this absorption is commonly 
seen in Type la SN spectra near maximum light, al¬ 
though it is also prominent in the spectra of Type lb and 
Ic core collapse SNe. The red color of the HFFMTom 
SED is qualitatively consistent with a redshift of z > 1 
- although this information was not used by SNID for 
the template matching. As we will see in Section [5j this 
spectral classification of SN HFFMTom is reinforced by 


the photometric information, which also supports classi¬ 
fication as a Type la SN at z~ 1.35. 

4.2. Spectral Fitting with UV Type la Templates 

To refine the redshift and phase constraints on 
HFFMTom, we next fit the spectrum with a custom spec¬ 
tral matching program that employs a library of Type la 
SN SEDs. This library is similar to the Type la spec¬ 
tral set used by SNID, but also includes more recent SNe 
with well-observed spectral time series that extend to 
rest-frame UV wavelengths (e.g. SN 201 lfe and 2014J). 
We first use an approach similar to the SNID algorithm: 
warping the pseudo-continuum of each template spec¬ 
trum by dividing out a 3 ld -order polynomial to match 
the observed SED of SN HFFMTom. This approach will 
find templates that have similar abundances and photo- 
spheric velocities. With the redshift fixed at z = 1.3457 
we find the best fit is a sp ectrum from the normal Type 
la SN 2014J ( |Foley et al.|2014| ), with a x 2 P er degree of 
freedom v of yfi jv — 121.8/89 shown in the top panel 
of Figure [2j The excess variance in this fit may be at¬ 
tributed to the intrinsic variation of Type la SN sp ectra , 
which is more prominent at UV wavelengths (e.g. Foley 
et al.|2008 Wang et al.|2012 ) and is not fully represented 
in the available template library. Allowing the redshift 
as an additional free parameter, we still find results that 
are consistent with the SNID fits and < 2<r from the host 
galaxy redshift: z = 1.31 ±0.02 and a phase of 0±3 rest- 
frame days. The best-fitting spectral template in th is 
case i s the normal Type la SN 2012cg (Amanullah et al. 
2015|) at= 1.31 with yfi/v = 98.5/88 (second panel of 
Figure [2]) . 

Next, we repeat the fitting, but without any warping 
of the templates to account for differences in the con¬ 
tinuum shape. In this iteration we only allow each tem¬ 
plate SED to be scaled in flux coherently at all wave¬ 
lengths, so the fits are more sensitive to the overall color 
of the SED P^I Fixing the redshift to z = 1.3457, we find 
the best match is fr om the normal Type la SN 201 lfe 
(Mazzali et al. |2014 ), though the fit is quite poor, with 
X 2 l v = 313.2/91 (third panel of Figure E]). When the 
redshift is allowed as a free parameter the HFFMTom 
SED is still matched best by a SN 201 lfe template, now 
at redshift z = 1.31 with \ 2 /i/ = 240.8/90 (fourth panel 
of Fi gure[2|). SN 2011fe had effectively n o dust reddening 
(e.g. Nugent et al.|2011 Li et al.|2011a), so the fact that 
SN 2011te provides the best un-warped template match 
is further evidence that SN HFFMTom suffers from very 
little dust extinction. 

Without the continuum warping, an alternative fit also 
arises: the fa st-declining (91bg-like) Type la SN 201 liv 


(Foley et al. 2012) at z = 0.98 ± 0.01. Formally, this 
match provides a slightly better fit to the unwarped 
HFFMTom spectrum (x 2 / v = 220.6/90, bottom panel of 
Figure |2J) , although the fit is notably poorer at ~ 8700A 
where the most significant absorption feature is found. 
Furthermore, a redshift z ~ 1 is at odds with the spec¬ 
troscopic redshift of the nearest galaxy (z = 1.3457), and 
we will see in the following section that the photometric 
data is also incompatible with a Type la SN at z ~ 1.0. 

39 Note that gravitational lensing does not affect the color of 
background sources at all, so the expected lensing magnification of 
HFF14Tom does not affect this analysis. 
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Setting aside the z ~ 1 solution, all other tem¬ 
plate matches provide a consistent redshift constraint of 
z = 1.31 ± 0.02, regardless of whether the templates are 
warped to match the SN HFF14Tom continuum shape. 
The inferred age from these fits is 0 ± 3 rest-frame days 
from peak brightness, which is also consistent with the 
observed light curve. Taken together with the host 
galaxy spectroscopic redshift, these fits suggest that SN 
HFF14Tom is a normal Type la SN at z = 1.3457 with 
an SED color close to SN 201 lfe, but spectral absorption 
features similar to SN 2014J. 


5. PHOTOMETRIC CLASSIFICATION 

Relative to other SNe at z > 1, the SN HFF14Tom 
light curve was unusually well sampled at rest-frame ul¬ 
traviolet wavelengths, due to the rapid cadence of the 
HFF imaging campaign. These ACS observations there¬ 
fore provide a tight constraint on the time of peak bright¬ 
ness and the evolution of the SN color. Supplemental 
observations with the WFC3-IR camera provided critical 
rest-frame optical photometry, enabling a measurement 
of the apparent luminosity distance through light curve 
fitting. 

As a check on the spectral classification of SN 
HFF14Tom (Section 4.1), we independently classified the 
SN using a Bayesian photometric classifier. We use the 
sncosmo software packagep*] to simulate SN light curves 
from 2 = 0.3 to 2.3 and evaluate the classification prob- 
ability using traditional Bayesian model selection (as in 


Jones et aLf201 3; Rodn ey et al.||2014| |Graur et al.||2014[ 


Rodney et al. 2015|). In this analysis we re present nor¬ 

mal T ype la SlNe with the SALT2 model (Guy et al. 
2010), and core collapse SNe with 42 discrete templates 
(26 Type II and 16 Type Ib/c) drawn from the tem¬ 
plate library o f the SuperNova Analysi s software pack¬ 
age (SNANA, Kessler et al. 2009b Ip*] Likelihoods are 
defined by comparing the observed fluxes to model pre¬ 
dictions in all passbands where the model is defined. In 
practice, this means we exclude the SN detections in the 
F435W and F606W bands, which are too blue for our 
models at z > 0.85. 

The core collapse SN models have free parameters for 
date of peak brightness (f p k)> amplitude, and redshift 
(z). Due to the expected impact of gravitational lensing 
magnification, we use a flat prior for the intrinsic lumi¬ 
nosity of all SN sub-classes. We also apply a flat prior 
for the SN redshift, which allows our photometric anal¬ 
ysis to provide an independent check on the host g alax y 
photo-z and spectroscopic redshift (Sections [3] and 14.21) . 

For Type la SNe, the SALT2 model has two additional 
parameters that control the shape (or) and color (c) of 
the light curve. We use conservative priors here, de¬ 
fined to encompass a range of Type la SN shapes and 
colors that is broader th an typically allowed in cosmo¬ 
logical analyses (see e.g., |K essler et al.||2009al [Sullivan 
et ah] |2011 ; Rest et al.|2014 FT For x\ the prior is a bifur¬ 
cated Gaussian distribution with mean x\ =0, disper¬ 
sion cr+ =0.9 and a~ 1 = —1.5. The bifurcated Gaussian 
prior for the color parameter c has c = 0.0, a~ = 0.08, 
and < 7 + = 0.54. The c parameter in SALT2 combines 
intrinsic SN color and extinction due to dust, so the 

40 http://sncosmo.github.io/ 

41 Throughout this work we use SNANA v!0_35g. 
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Figure 3. Maximum likelihood model for each SN sub-class, 
derived from Bayesian model selection using the photometric data 
alone. Grey points show the observed SN HFF14Tom photometry 
with error bars, though these are typically smaller than the size of 
the marker. The Type la model (orange solid line) is drawn from 
the SALT2 template at z = 1.35. The best match from all Type 
Ib/Ic models is based on the Type Ic SN SDSS-14475 at z = 0.695 
(green dashed line). For the Type II class, the best match is from 
the Type II-L SN 2007pg at z = 1.8 (blue dash-dot line). The 
Type la model is by far the best match, and the only one that 
is consistent with both the spectroscopic redshift of the probable 
host galaxy and the spectroscopic redshift from the SN spectrum. 
The date of the HST spectral observations is marked with a thin 
grey vertical line. 


large red tail of this distribution allows for the possi¬ 
bility of several magnitudes of dust extinction along the 
HFF14Tom line of sight. 

We also assign a class prior for each of the three pri¬ 
mary SN sub-classes (Type la, Ib/c, and II), using a fixed 
rel ative fraction f o r each sub- class as determine d at z = 0 
by Smartt et al. (2009) and Li et al. (2011b). A more 
rigorous classification would extrapolate these local SN 
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class fractions to higher redshift using models or mea¬ 
surements of the volumetric SN rate. For simplicity, we 
do not vary the class priors with redshift, and in practice 
these priors do not have any significant impact on the 
resulting classification. 

The final photometric classification probability for SN 
HFF14Tom is p(Ia|D) = 1.0, with the classification prob¬ 
ability from all core collapse SN sub-classes totaling less 
than 10 -32 . Although this Bayesian classification utilizes 
the full posterior probability distribution, for illustration 
we highlight in Figure [3] a single best-fit model for each 
sub-class. This demonstrates how the core collapse SN 
models fail to adequately match the observed photome¬ 
try. In particular, only the Type la model can simultane¬ 
ously provide an acceptable fit to the well-sampled rising 
light curve in F814W and the F814W-F160W color near 
peak. 

The marginal posterior distribution in redshift for the 
Type la model is sharply peaked at 0 = 1.35 ± 0.02, 
which is fully consistent with the redshift of the pre¬ 
sumed host galaxy {z = 1.3457) as well as the redshift 
of z = 1.31 ± 0.02 derived from the SN spectrum in Sec¬ 
tion |4.2| The time of peak brightness is also tightly con¬ 
strained at f p k = 56816.3 ± 0.3, which means the HST 
grism observations were collected within 2 rest-frame 
days of the epoch of peak brightness - also consistent 
with our spectroscopic analysis. 


6. DISTANCE MODULUS AND MAGNIFICATION 

With the type and redshift securely defined as a nor¬ 
mal Type la SN at z = 1.3457, we now turn to fitting the 
light curve with Typ e la templates to measure the dis¬ 
tance modulus ()6.1) a nd t hen derive the gravitational 
lensing magnification (16.2). In this process, we would 
like to avoid introducing systematic uncertainties inher¬ 
ent to an y ass umed cosmological model. To that end, 
in section |6.2| we follow P14 and define the magnifica¬ 
tion by comparing the measured distance modulus of SN 
HFF14Tom against an average distance modulus derived 
from a “control sample” of unlensed Type la SNe at sim¬ 
ilar redshift. This allows us to make only the minimal as¬ 
sumption that the redshift-distance relationship for Type 
la SN is smooth and approximately linear over a small 
redshift span, which should be true for any plausible cos¬ 
mological model. 


6 .1. Light Curve Fitting 

As in m we derive a distance modulus for SN 
HFF14Tom and all SNe in our control sample using two 
independent light curve fitters: the S ALT2 model de¬ 
scribed above and the MLCS2k2 model (Jha et al.|2007). 
With both fitters we find light curve shape and color pa¬ 
rameters for HFF14Tom that are fully consistent with 
a normal Type la SN. For SALT2, with the redshift 
fixed at z = 1.3457, we find a light curve shape pa¬ 
rameter of x\ = 0.135 ± 0.199 and a color parameter 
of c = —0.127 ± 0.025, yielding a \ 2 value of 45.0 for 
36 degrees of freedom, v. With the MLCS2k2 fitter the 
best-fit shape parameter is A = —0.082 ± 0.070 and the 
color term is Ay = 0.011 ±0.025, giving % 2 jv = 22.8/36. 

The MLCS2k2 fitter returns a distance modulufF 2 ! 
dmMLCS 2 k 2 directly, as it is defined to be one of the free 


42 We use ’dm’ to indicate the distance modulus to avoid con- 


Table 3 

HFF14Tom Measured Distance Modulus and 
Magnification at z = 1.3457 



Distance Modulus 

Measured 

Fitter 

HFF14Tom 

Control 

Magnification 

MLCS2k2 

44.205 ±0.12 

44.97 ±0.06 

2.03 ±0.29 

SALT2 

44.177 ± 0.18 

44.92 ± 0.06 

1.99 ±0.38 


parameters in the model. To derive dm from the SALT2 
fit, we use 


dmsALT 2 = rn* B — M + a(s - 1) — f3C. 


( 1 ) 


Here the parameters for light curve shape s and colo r C 
corre spond to the SiFTO light curve fit ter (Conley et al . 
20081, so we first use the formulae from [Guy et al. ifSOTOf 
to convert from SALT2 ( x\ and c) into the equivalent 
SiFTO parameters. We also add an offset of 0.27 mag 
to the value of m* B returned by SNANA, in order to 
match t he arbitrary norma lizati on of the SALT2 fitte r 
used by |Guy et al.| d2010| and [Sullivan et al.| (|2011|). 
This conversion from SALT2 to SiFTO is necessary, as 
it allo ws us to a dopt va lues for the constants M, a , and 
/3 from Sullivan et al. (2011), which have been calibrated 


using 472 SJNe from the SNLS3 sample (Conley et al. 


2011): M = -19.12 ± 0.03, a = 1.367 ± 0 086, and ft = 
TT79± 0.101. 

The SNANA version of the MLCS2k2 fitter returns 
a value for the distance modulus (dmMLCS 2 k 2 ) that has 
an arbitrary zero point offset relative to the SALT2 dis¬ 
tances (dmsALT 2 )- To put the two distances onto the 
same reference frame we add a zeropoin t cor rection of 
0.20 mag to the MLCS2k2 distances as in |P14[ This cor¬ 
rection was derived by applying both fit ters to a sam ple 
of Type la SNe from the S PSS survey (Holtzman et al. 


2008 


Kessler et al. 2009a), with the extinction law Ky 


fixed at 1.9. 

The total uncertainty in the distance modulus is 


^tot — 


+ < J tnf 


( 2 ) 


The cj s tat term is the statistical uncertainty, which en¬ 
capsulates uncertainties from the data and the model, 
and (Ti n t accounts for the remaining unmodeled scatter. 
This latter term is derived by finding the amount of addi¬ 
tional distance modulus scatter that needs to be added to 
a Type la SN population to get x 2 P er degree of freedom 
equal to 1 for a fiducial cosmological model fit to the SN 
la Hubble diagram (distance modulus vs redshift). Thus, 
(Tint is designed to account for any unknown sources of 
scatter in the Type la SN population, including uniden¬ 
tified errors in the data analysis as well as the natural 
scatter in intrinsic Type la SN luminosities. The value 
of (Ti nt may be expected to vary as a function of redshift 
and also from survey to survey. 

Recent work has highlighted the inadequacy of this 
simplistic approach for h andling intrinsic scatter in the 
Type la SN population (Marriner et al. 2011 Kessleij 


|et al.|[2013| |Mosher et al.||2014| |Scolnic et al.||2() 14b| |Be* 


fusion, reserving the symbol ^ to refer to the lensing magnifi¬ 

catio n" This ’am’ is a standard distance modulus, defined as 
dm = 51og 10 cLl + 25, where d l is the luminosity distance in Mpc. 
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Figure 4. Type la light curve fits to SN HFF14Tom using the MLCS2k2 (top row) and SALT2 (bottom row) fitters. The model redshifts 
are set to 2 : = 1.3457 as determined from combined spectroscopic and photometric constraints. Solid lines denote the best-fit model and 
shaded lines show the range allowed by 1-cr uncertainties on the model parameters. Observed fluxes are shown as diamonds, scaled to an 
AB magnitude zero point of 25. Error bars are plotted, but most are commensurate with the size of the points. The left-most panel includes 
observations in the F435W and F606W filters, although these were not used for the fit, as they are bluer than the minimum wavelength 
for the model. The lower axis marks time in observer frame days, while the top axis shows the time in the rest frame relative to the epoch 
of peak brightness. 


toule et aL[|20141, but a full consideration of those alter¬ 


native approaches is beyond the scope of this work. We 
adopt the simple approach of using a single empirically 
defined value for (Tint, reflecting principally an intrinsic 
scatter in Type la SN intrinsic luminosity that is fixed 


0.08 mag (Jha et al. 

2007; Conley et al. 

2011) to 0.15 

mag ( 

Kessler et al. 2009a 

Suzuki et al. 201 

12). We adopt 

a value of a iut = 0.08 mag, as derived 

by |Jha et al. 


(2007) and Conley et al. (2011). Although this is on 
the low end of the range reported in the literature, this 
value is the most appropriate to apply to our analysis 
for two reasons. First and foremost, for the SALT2 fit¬ 
ter we are using light curve fit parameters (a,/3) with 
associated uncertainties that have been derived from the 


joint_analysis of Conley et al. (2011) and Sullivan et al. 
(2011). Similarly, the implementation of MLCS2k2 that 
we have used is based on the uncertainty model derived 
from model training in Jha et al. (2007). Inflating Omt 
beyond 0.08 mag would therefore be equivalent to driv¬ 
ing the reduced y 2 of the Type la SN Hubble diagr am to 
< 1. Second, th e value of 0.08 mag determined in |Con- 


ley et al. 


Ml i 

I (2007 


is specific to the HST SN sample, from 


Kiess et al. 1 20071 and Suzuki et al. (2012), which is the 
SIN subset that is most similar to HI F 14Tom in terms of 
redshift and data analysis. Larger values for <7i nt are typ¬ 
ically associated with SN samples at significantly lower 
redshifts that have less homogeneous data collection and 
analysis. 

Final values for the SN distance modulus are shown in 
Table[3j adopting the spectroscopic redshift (z = 1.3457). 


In addition to modifying the distance modulus uncer¬ 
tainty for SN HFF14Tom, we also add dint = 0.08 mag 
in quadrature to the uncertainty for every SN in the con¬ 
trol sample. Note however that the uncertainty on the 
control sample value at z = 1.3457 is only 0.06 mag, 
smaller than the intrinsic dispersion of any single SN 
because it reflects our measurement error on the mean 
distance modulus of the population. The distance mod¬ 
uli derived from the SALT2 and MLCS2k2 light curve 
fitters are fully consistent within the uncertainties. 


6.1.1. Host Galaxy Mass Correction 

It is now an accepted practice in cosmological analyses 
using Type la SN to apply a correction to the luminosity 
of each SN based on the stellar mass of its host galaxy. 
Typically this is described as a simple bifurcation of the 
SN population: SNe that appear in more massive hosts 
are observed to be ~0.08 mag brighter (after corrections 
for lig ht curve shape and color) than SNe in low-mass 
hosts ( Kelly et al.|2010 Sullivan et al.poTO I. The divid¬ 
ing line tor this purely empirical “mass step” correction 
is generally set around 10 lo M g . Although this t h resh old 
value is somewhat arbitrary (see e.g. Betoule et al.|2014), 
it happens to be very close to the -SIN HFKl lTom host 
galaxy mass of 10 9 8 Mq (Section EJ). 

The physical mechanism that arives the mass step is 
not yet understood, but may be related to the metallicity 
or age of the SN progenitor systems. In either case, the 
significance of this effect should decrease with redshift, 
as metal-rich passi ve galaxies become much less common 
at z > 1 (see e.g. Rigault et al. 2013 Childress et al. 
2014|. Indeed, when the size of the mass step correction 
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is allowed to vary with redshift, there is no significant 
evidence that a non-zero mass step is required at z > 1 


( Suzu ki et al.|2012 [Shafer & Huterer|2014| |Betoule et al. 
1 the abi 


2014| ). Given 


sence of a clear physical model and 
the lack of empirical support for a high -z mass step, we 
do not apply any correction to SN HFFMTom or the 
control sample. 

If we were to apply the correction to this sample using 
the standard approach, then the SN HFFMTom distance 
modulus would not be adjusted, because its host galaxy 
mass (Section [3]) is just below 10 10 M 0 . Note that there 
is some circularity here, as we used the SN magnifica¬ 
tion in deriving the host galaxy mass, and now use the 
host mass to inform the magnification. However, regard¬ 
less of whether the HFFMTom host mass falls above the 
mass threshold or below, this correction is not substan¬ 
tial enough to account for the observed magnification 
tension. The appropriate correction to be applied to a 
SN above the mass threshold in this redshift range would 
be <0.04 mag ( Rigault et al.|2013 Childress et al.|2014 ), 
which results in a change to the interred magnification of 
HFFMTom that is much less than the observed discrep¬ 
ancy of ^0.24 mag. 



1.2 1.3 1.4 1.5 

Redshift 


6 .2. Control Sample Comparison 

The unlensed sample comprises 22 spectroscopically 
confirmed Type la SNe in the range 1.1 < z < 1.6 


z SN 
7 fror 
DELS 

searcl 
n SCP 

43 

vt 

(Strolger et al. 2004 

Riess et al. 

2007). 

i Suzuki et al.j 2012), and 4 from CAN- 

(Rodney et al. 2012 2014 

). Using the SALT2 


and MECS2k2 litters as described aDOve, we get distance 
modulus measures for every object in this control sam¬ 
ple. We then fit a linear relationship for distance modu¬ 
lus vs. redshift, and derive a prediction for the distance 
modulus of a normal Type la SN at the redshift of SN 
HFFMTom (Figure [5|. This predicted value is given in 
Table [3] under the “Control” column. The difference be¬ 
tween the observed distance of SN HFFMTom and this 
control sample value is attributed to the magnification 
from gravitational lensing: 


dm contro i - dm H FFi4Tom = 2.5 log 10 p. (3) 

The inferred magnifications for the two different fitters 
are reported in the final column of Table [3] Although we 
have chosen to use a cosmology-independent approach 
to determine the magnifications, we note that these re¬ 
sults are fully consistent with the values that would be 
determined by comparing the observed SN HFFMTom 
distance modulus against the predicted value from a flat 
ACDM cosm ology . Adopting the cosmological parame- 


ters used in P14 and derived in Sullivan et al. (2011) 
(Hq=71.6, fl rn =().27. Ha=0.73) would give a predicted 
distance modulus dmACDM = 44.90 and a measured 


logical parameters (e.g., 

Betoule et al. 2014 

Planck Col- 

laboration et al. 2015 

) would not significant 

y change the 


43 GOODS: the Great Observatories Origins Deep Survey, 
PLGiavalisco, HST-PID:9425,9583 

44 SCP: the Supernova Cosmology Project, PLPerlmutter 

45 CANDELS: the Cosmic Assembly Near-infrared Deep Extra- 
galactic Legacy Survey, PLFaber & Ferguson 


Figure 5. Measurement of the lensing magnification from com¬ 
parison of the HFF14Tom distance modulus to a sample of un¬ 
lensed field SN. SNe from the GOODS program are plotted as 
squares, those from the SCP survey are shown as circles, and the 
CANDELS objects are triangles. The distance modulus for each 
SN is derived from light curve fits using the MLCS2k2 fitter (top 
panel) and the SALT2 fitter (bottom panel). Grey lines with shad¬ 
ing show a linear fit to the unlensed control sample is shown, an¬ 
chored at the redshift of SN HFF14Tom, with fit parameters given 
at the top. The derived magnification is reported at the bottom of 
each panel. 


7. COMPARISON TO MODEL PREDICTIONS 

Before the Frontier Fields observations began, the 
Space Telescope Science Institute (STScI) issued a call 
for lens modeling teams to generate mass models of all 
6 Frontier Field clusters, using a shared collection of all 
imaging and spectroscopic data available at the time. In 
response to this opportunity, five teams generated seven 
models for Abell 2744. These models necessarily relied 
on pre-HFF data, and were required to be complete be¬ 
fore the HFF program began, in order to enable the esti¬ 
mation of magnifications for any new lensed background 
sources revealed by the HFF imaging. An interactive 
web tool was created by D. Coe and hosted at STScI, to 
extract magnification estimates and uncertainties from 
each model for any given redshift and position. In this 
work we also consider 8 additional models that were cre¬ 
ated later. Some of these are updates of the original 
models produced in response to the lens modeling call, 
and some of them include new multiply-imaged galaxies 
discovered in the HFF imaging as well as new redshifts 
for lensed background galaxies. Magnifications derived 
from several of these later models are also available on 
the interactive web tool hosted by STScI. The full list 
of models and details on their construction are given in 
Table [4j 

Table[5]gives each model’s predicted magnification and 
uncertainty for a source at z = 1.3457 and at the posi¬ 
tion of HFFMTom. For thirteen of these models the SN 
magnification qualifies as a true “blind test”, as they were 
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Table 4 

Tested lens models for Abell 2744. 


Model N sys a Ni m b N spe c c N p hot d References Description 


Bradac(vl) e 

16 

56 

2 

14 

HFF f 

Bradac et al. 2009 

SWUnited g : Free-form, strong+weak-lensing 







based model.Errors from bootstrap resampling 
only weak-lensing constraints. 

CATS(vl) 

17 

60 

2 

14 

HFF; 

Richard et al. 2014 

CATS b team implementation of LENSTOOL 1 
parametric strong-lensing based model. 

CATS(vl.l)J 

17 

60 

2 

14 

Richard et al. |2014 

CATS team LENSTOOL parametric model with 
both strong and weak lensing constraints. 

Merten 6 

16 

56 

2 

14 

HFF; 

Merten et al. 2011 

SaWLENS, k Grid-based free-form strong-(-weak 
lensing based model using adaptive mesh 
refinement. 

Sharon(vl) 

17 

60 

2 

14 

HFF; 


LENSTOOL parametric strong-lensing based 
model 

Sharon(v2) 

15 

47 

3 

11 

Johnson et al.|2014| 

LENSTOOL parametric strong-lensing based 
model. Includes cosmological parameter varia¬ 
tions in uncertainty estimates. 

Zitrin-LTM 

10 

44 

2 

0 

HFF; 

Zitrin et al.|2009 

Parametric strong-lensing model, adopts the 
Light-Traces-Mass assumption for both the lumi¬ 
nous and dark matter. 

Zitrin-NFW 

10 

44 

2 

0 

HFF; 

Zitrin et al.|2013| 

Parametric strong-lensing model using PIEMD 
profiles for galaxies and NFW profiles for dark 
matter halos. 

Williams 

10 

40 

2 

8 

HFF 


GRALE 1 : Free-form strong-lensing model using a 
genetic algorithm. 


Post-HFF models : include data from the HFF program 


CATS(v2) 

50 

151 

4 

1 |Jauzac et al.|2015| 

Updated version of the CATSvl model, adds 33 
new multiply-imaged galaxies, for a total of 159 
individual lensed images. 

Diego m 

15 

48 

4 

11 Diego] 2014 

WSLAP+ n : Free-form strong-lensing model using a 
grid-based method, supplemented by deflections 
fixed to cluster member galaxies. 

GLAFIC 

24 

67 

3 

12 Ishigaki et al. ]2015| 

Parametric strong-lensing model using vl.O of the 
GLAFIC code. 0 

Lam(vl) 

21 

65 

4 

17 Lam et al. |2014| 

Alternative implementation of the WSLAP+ model, 
using a different set of strong-lensing constraints 
and redshifts. 


Unblind models : generated after the SN magnification was known 


Bradac(v2) e 

25 

72 

7 

18 

Wang et al. |2015| 

Updated version of the SWUnited model with new 
strong-lensing constraints from HFF imaging and 
GLASS spectra. Errors derived via bootstrap re¬ 
sampling only strong-lensing constraints. 

Lam(v2) 

10 

32 

5 

5 

Lam et al. |2014| 

Updated version of the WSLAP+ model, using 
more selective strong-lensing constraints and 
GALFIT P models for galaxy mass. 

CATS(v2.1) 

55 

154 

8 

1 

pauzac et al.|2015| 

Updated version of the CATSv2 model, adopting 
the spec-z constraints used for Bradac(v2). 

CATS(v2.2) 

25 

72 

8 

1 

pauzac et al.|2015| 

Updated version of the CATSv2 model, adopting 
the spec-z constraints and multiple-image defini¬ 
tions used for Bradac(v2). 


a Number of multiply imaged systems used as strong-lensing constraints. 
b Total number of multiple images used. 

c Number of multiply imaged systems with spectroscopic redshifts. 
d Number of multiply imaged systems with photometric redshifts. 

e Uses weak lensing constraints from HST- and ground-based imaging, as described in [Merten et al.| ( |201l| . 

f Lens models with the reference code “HF F” were produced a s part of the Hubble Frontie r Fields lens modeli ng program, using 
arcs identified in HST archival imaging from [M erten et al.l2011 | spectroscopi c redshifts fromjRichard et al. 2014, and ground-based 
imaging from |Cypriano et al7||2QQ-4| |(Jkabe &; Um etsu 200H| fOkabe et al.||2UlUa|b| Details on the model construction and an 
interactive model magnification web interface are a vailable at httpT77archive.stsci.edu/prepds/frontier/lensmodels/ 


g SWunited : Strong and Weak lensing United; |Bradac et al.|2UU5 


CATS : Clusters As TelescopeS lens modeling team. Pi’s: J.-P. Kneib & P. Natarajan 
1 LENSTOOL : |Julio et al.|2007| http://projects.lam.fr/repos/lenstool/wiki 
j Uses weak lensing constraints from HST imaging, as described in |Richard et al7| (|2014[) 
k SaWLENS : |Merten et al.|2009[ St rong and Weak LENSing analysis code. http://www. 

1 GRALE : GRAvitational LEnsing;|Liesenborgs et al.|2006||2007||T 


Mohammed et al. 2014] 
Abell 2744 model available at http://www.ifca.unican.es/users/jdiego/LensExplorer 


ulianmerten.net/codes.html 


No uncertainty estimates were avail¬ 
able for th e Diego i mplement ation of the WSLAP+ model, so we adopt the uncertainties from the closely related Lam model. 
n WSLAP+ : iSendra et al.|2014[ Weak and Strong Lensing Analysis Package plus member galaxies (Note: no weak-lensing constraints 


used for Abell 2744) 


GLAFIC : Oguri 2010; http://www.slac.stanford.e du/~oguri/glafic / 
GALFIT : Two-dimensional galaxy fitting algorithm ( |Peng et al |2002| ) 
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Table 5 

Lens model predictions for SN HFF14Tom 
magnification 


Model 3, 

Best b 

Median c 

68% Conf. Range d 

Bradac(vl) 

3.19 

2.48 

2.31-2.66 

CATS(vl) 

2.28 

2.29 

2.25-2.34 

CATS(vl.l) 


2.62 

2.44-2.80 

Merten 

2.33 

2.24 

2.04-2.92 

Sharon(vl) 

2.56 

2.60 

2.44-2.78 

Sharon(v2) 

2.74 

2.59 

2.42-2.85 

Zitrin-LTM 

2.67 

2.99 

2.61-3.77 

Zitrin-NFW 

2.09 

2.29 

2.07-2.52 

Williams 

2.70 

2.81 

1.65-5.54 

CATS(v2) 


3.42 

3.27-3.58 

Diego 


1.80 

1.44-2.16 

GLAFIC 

2.34 

2.29 

2.19-2.37 

Lam(vl) 


2.79 

2.42-3.16 

Bradac(v2)* 

2.23 

2.26 

2.30-2.23 

Lam(v2)* 

1.86 

1.91 

1.54-2.28 

CATS(v2.1)* 


3.06 

2.92-3.19 

CATS(v2.2)* 


3.07 

2.94-3.20 


3 Models above the line are from the pre-HFF set, and 
those below incorporate HFF data. The final four, 
marked with asterisks, were not formally part of the 
blind test as they include modifications made after the 
measured magnification of the SN was known. 
b The magnification returned for the optimal version 
of each model, as independently defined by each lens 
modeling team. 

c Median magnification from 100-600 Monte Carlo re¬ 
alizations of the model. 

d Confidence ranges about the median, enclosing 68% 
of the realized values. 

completed before the SN HFF14Tom magnification mea¬ 
surement was known. The final four models (Bradac-v2, 
Lam-v2, CATSv2.1 and 2.2) are technically not blind, al¬ 
though none of these modelers used the SN magnification 
as a constraint, and the modelers did not have access to 
the final SN magnification value when constructing their 
model. 

These models represent a broad sampling of the tech¬ 
niques and assumptions that can be applied to the mod¬ 
eling of mass distributions in galaxy clusters. We discuss 
two primary modeling choices here, but a more rigorous 
comparison of these diverse lens modeling techniques is 
beyond the scope of this work. For a complete discussion 
of each model’s methodology the reader is directed to the 
listed references. 


(1) Parametric vs. Free-form: Ten of these models 
are so-called parametria ^ q models, and the remaining 
seven are free-form models. Broadly speaking, the para¬ 
metric models use parameterized density distributions 
to describe the arrangement of mass within the cluster. 
Therefore, parametric models rely (to varying degrees) 
on the assumption that the cluster’s dark matter can be 


descri bed by analytic forms such as NFW halos (Navarro 


et al. 1997|), or p seudo isothermal elliptica l mass distri¬ 
butions (FlEMD Kassiola & Kovner||1993). 

The seven free-iorm models divide the cluster field into 
a grid, generally using a multi-scale grid to get better 


46 Although this nomenclature is becoming standard in the liter¬ 
ature, it is somewhat misleading, as the pixels or grid cells in free¬ 
form models are effectively parameters as well. Perhaps “simply- 
parameterized” would be more accurate, though we adopt the more 
common usage here. 


sampling in regions with a higher density of information 
(e.g. density of multiple images). Each grid cell is as¬ 
signed a mass or a potential, and then the mass values 
are iteratively refined to match the observed lensing con¬ 
straints. In some cases an adaptive grid is used so that 
the grid spacing itself can also be modified as the model 
is iterated (e.g. Liesenbo rgs et al. 2006 Merten et al. 
2009 |Bradac et al.||2(l(j9| ~ 


"As usual, there is a tradeoff between a model’s flex¬ 
ibility, the strength of the model assumptions, and the 
resulting uncertainties. In a probabilistic framework, the 
posterior distribution function of the desired quantities 
depends on all the priors, including model assumptions 
like parametrization. In general, free-form methods tend 
to be more flexibile than simply parameterized models, 
and thus tend to result in larger error bars. In brief, if the 
free form methods are too flexible, then they will result 
in overestimated error bars. Conversely, if the simply 
parameterized models are too inflexible, they will result 
in underestimated error bars. 

(2) Strong vs. Strong + Weak: Thirteen of the mass 
models evaluated here are constrained only using strong- 
lensing features such as multiply-imaged background 
galaxies and highly magnified arcs. There are now 


150 know n lensed images behind Abell 2744 (Jauzac 
et al. 2015), but SN HFF14Tom is located several arc- 
seconds outside the core region of the cluster where 
these multiply-imaged galaxies are found. Therefore 
these “strong-lensing only” models must necessarily rely 
on extrapolations to provide a prediction for the SN 
HFF14Tom magnification. 

The other four models in our comparison set also use 
the strong lensing constraints from the core region, but 
additionaly use weak lensing measurements to provide 
additional constraints on the mass distribution farther 
from the cluster core. Notes on the derivation of these 
weak-lensing constraints are given in the final column 
of Table [3] The signal from weak lensing relies on a 
large sample of background galaxies, so this constraint 
operates principally at separations more than V from 
the cluster core. At a projected separation of ~ 40", SN 
HFF14Tom falls in between the strong- and weak-lensing 
regimes. 

In Figure [6] the model predictions are plotted along¬ 
side the observed magnification of SN HFF14Tom, de¬ 
rived in Section [6] To first order, this comparison 
shows that these 17 models are largely consistent with 
each other and with the observed magnification of SN 
HFF14Tom. The “naive mean” of the full set of models 
is /r pre = 2.6 ± 0.4. This is an unweighted mean (i.e. 
we ignore all quoted uncertainties) derived by naively 
treating each as an independent prediction for the mag¬ 
nification (this is clearly incorrect, as several models are 
represented multiple times as different versions). The 
naive mean is separated from the observed SN magnifi¬ 
cation by Sgi/p = 28%, which is approximately a 1.5cr 
difference. 

This is approximate l y con sistent with the results of 


P14 and Nordin et al. (2014), where model predictions 
were found to be in reasonably good agreement with a set 
of 3 lensed SNe from the CLASH program. The general 
agreement between the model predictions and the SN 
measurement is especially encouraging for these Abell 
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Figure 6. Comparison of the observed lensing magnification to 
predictions from lens models. The vertical blue line shows the 
constraints from SN HFF14Tom derived in Section [6] using the 
MLCS2k2 fitter, with a shaded region marking the total uncer¬ 
tainty. Markers with horizontal error bars show the median mag¬ 
nification and 68% confidence region from each of the 17 lensing 
models. Circles indicate models that use only strong-lensing con¬ 
straints, while diamonds denote those that also incorporate weak- 
lensing measurements. Models using a “free-form” approach are 
shown as open markers, while those in the “parametric” family 
are given filled markers. The top half, with points in black, shows 
the nine models that were constructed using only data available 
before the start of the Frontier Fields observations. The lower 
eight models in green used additional input constraints, includ¬ 
ing new multiply imaged systems and redshifts. The final four 
points, with square orange outlines, are the “unblind” models that 
were generated after the magnification of the SN was known. The 
black dashed line marks the unweighted mean for all 17 models, at 

[i = 2 . 6 . 


2744 models. This is a merging cluster with a complex 
mass distribution, and the SN is located outside of the 
strong-lensing region where the models are most tightly 
constrained. 

However, beyond this first-order agreement, there is 
a small systematic bias apparent. All but two of the 
lens models return median magnifications that are higher 
than the observed value, and six of the models are dis¬ 
crepant by more than 1.5er. These six discrepant models 
are all biased to higher magnifications. They are found in 
both the pre-HFF and post-HFF models, in the paramet¬ 
ric and free-form families, and among the strong-lensing- 


only and the strong+weak subsets. It is important to 
emphasize that SN HFF14Tom only samples a single line 
of sight through the cluster, and this bias to higher mag¬ 
nifications is minor. Nevertheless, a systematic shift of 
this nature is surprising, given the wide range of mod¬ 
eling strategies, input data, and physical assumptions 
represented by this set of models. In the following sub¬ 
sections we examine possible explanations for this small 
but nearly universal bias. We first consider whether a 
misinterpretation of the data on the SN itself can account 
for the observed systematic bias, and then examine the 
lens models. 


7.1. Possible Errors in Supernova Analysis 
7.1.1. Redshift Error 

If the redshift of the SN derived in Section [3] were in¬ 
correct, then one would derive a different value for the 
magnification, both from the SN measurement and the 
lens model predictions. Conceivably, this could resolve 
the tension between the measurement and the models. 
It is often the case in SN surveys that redshifts are as¬ 
signed based on a host galaxy association, typically in¬ 
ferred from the projected separation between the SN and 
nearby galaxies. In this case the redshift is strongly sup¬ 
ported by evidence from the SN itself: we find a con¬ 
sistent redshift from both the SN spectrum (Section 4.21 
and the light curve Section[5j which are both within lcr of 
the spectroscopic redshift for the nearest detected galaxy: 
z = 1.3457. This appears to be a solid and self-consistent 
picture, so the evidence strongly disfavors any redshift 
that is significantly different from z = 1.35. 

We have adopted the most precise redshift of 2 = 
1.3457 from the host galaxy as our baseline for the mag¬ 
nification comparison. If instead we adopt the spectro¬ 
scopic redshift from the SN itself (z = 1.31; Section 4.21 
then we find no significant change in the inferred magni¬ 
fications or in the suggestion of a small systematic bias. 

7.1.2. Foreground Dust and SN Color 

All SN sight-lines must intersect some amount of fore¬ 
ground dust from the immediate circumstellar environ¬ 
ment, the host galaxy, and the intergalactic medium 
(IGM). In the case of SN HFF14Tom one might posit 
some dust extinction from the intra-cluster medium 
(ICM) of Abell 2744, although measurements of rich clus¬ 
ters s uggest that the ICM has only a negligible dust con - 
tent (Maoz 1995 Stickel et al. 12002 Bai et al. 2007). 
When fitting the HFFI4Tom light curve we account tor 
dust by including corrections that modify the inferred lu¬ 
minosity distance based on the SN color. If after applying 
these dust corrections we were still underestimating the 
effect of dust along this sight-line, then the SN would 
appear dimmer than it really is, the inferred distance 
modulus would be higher, and the measured magnifica¬ 
tion would be biased to an artifically low value. Thus, an 
underestimation of dust would be in the right direction 
to match disc repa ncy we observe. 

In Section lfi.ll we found that SN HFF14Tom is on 
the blue end of the normal range of Type la SN colors. 
With the SALT2 fitter we measured a color parameter 
c = —0.127 ± 0.025, and with MLCS2k2 we found the 
host galaxy dust extinction to be Ay = 0.011 ± 0.025 
magnitudes. These colors are tightly constrained, as we 
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are fitting to photometry that covers a rest-frame wave¬ 
length range from ~ 3500 — 7000A and extends to ~30 
days past maximum brightness. This leaves little room 
for the luminosity measurement to be biased by dust, 
as the dimming of HFF14Tom would also necessarily be 
accompanied by some degree of reddening. 

Nevertheless, one might suppose that a bias could be 
introduced if we have adopted incorrect values for the 
color correction parameter (3 or extinction law Ry in the 
SALT2 and MLCS2k2 fits, respectively. The appropriate 
value to use for this color correction and how it affects 
inferences about the intrinsic scatter in Type la SN lu¬ 
minosities is a compl ex questi o n that is be yond the scop e 


of this work (see e.g. iMarr i ner et al.|2011 Chota rd et al. 
2011 Kessler et al. |2013[ ScofmceT^T 2014bp . How- 


ever, we can already rule this out as a solution for the 
magnification discrepancy. Our error on the HFF14Tom 
distance modulus already includes an uncertainty in the 
extinction law and a related error to account for the in¬ 
trinsic luminosity scatter. These are well vetted param¬ 
eters, bas ed on observations of ~ 500 SNe extending to 
z ~ 1.5 (Sullivan et al. 2011). Furthermore, there is 
no reason to propose that SN HFFMTom is uniquely af¬ 
fected by a peculiar type of dust. Thus, any change in the 
color correction applied to SN HFF14Tom would require 
the same adjustment to be applied to the unlensed SNe 
at similar redshift that make up our comparison sample, 
largely negating the effect on the inferred magnification. 

The traditional color corrections as formalized in SN 
light curve fitters are designed to account for a dust com¬ 
ponent that lies in the rest frame of the SN. The inferred 
luminosity of a SN can also be affected by the presence of 


different reddening law I 

Menard et al. 

2010 a) 

the magnitude of such a bias 

is insuij 

icient 

for the observed discrepancy, 

Menard et al. ( 

timate the opacity ot the universe as ( Ay) ^ 


0.03 mag 

up to z = 0.5. While this can have a measurable impact 
on precise cosmological constraints, it is far less than the 
0.23 mag difference between the observed magnification 
of HFF14Tom and the mean of the model predictions. 

Although the very blue color of SN HFFMTom is help¬ 
ful to rule out dust as a possible explanation for the ob¬ 
served small systematic bias, it is possible that this very 
blue color is itsel f leading to a bias in the distance mea- 

(2014a) measured a small bias 


surement. 


Scolnic et al. 


in the SAL'i’2 titter for Type la SNe that have a color 
parameter derived from the light curve fit c < —0.1 (see, 
e.g., their Figure 10 and Section 5.1). For HFFMTom 
we have measured c = —0.13, which would correspond 
to a bias of roughly +0.05 mag in the SALT2 distance 
measurement. If we were to apply a —0.05 mag cor¬ 
rection to the distance modulus from the SALT2 fit, 
then this would increase the inferred magnification to 
MSALT 2 = 2.08 ± 0.36. This would slightly reduce the 
tension between the SALT2 measurement and the model 
predictions, though it would not be enough to completely 
alleviate it. 


7.1.3. Misclassification 

Is it possible that HFFMTom is an example of a pe¬ 
culiar stellar explosion that does not follow the relation¬ 
ship between light curve shape and l umin osity observed 
in normal Type la SNe? In Sections |4.1| and [5] we clas- 
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Figure 7. Comparing HFF14Tom observations to the color 
curves of SN 2006bt and a normal Type la SN model. Grey points 
show the observed colors of SN HFFMTom, k-corrected to rest- 
frame UBVr bands. Solid orange lines show the observed color 
curves for the peculiar Type la SN 2006bt. Dashed blue lines show 
the color curves for a normal Type la SN, derived from the SALT2 
model with a color parameter c = —0.13, set to match the best-fit 
color for SN HFFMTom. 


sified SN HFFMTom as a normal Type la SN based on 
both the spectroscopic and photometric evidence. This 
rules out the possibility that HFFMTom belongs to a 
different class of n ormal SN explosions (Type lb, Ic, or 
II). In Section 6.1 we fit the HFFMTom light curve to 
determine a luminosity distance and found that the light 
curve shape and color are consistent with a Type la SN of 
average light curve width, with minimal dust extinction. 
This excludes the possibility that HFFMTom is one of 
the sub-class of faint and fast-declining Type la SNe like 
SN 1991bg. 

One remaining possibility is that HFFMTom could be 
part of a rare sub-category of peculiar Type la SNe that 
masquerade as their n ormal cousins, ep itomized by the 


prototype SN 2006bt (Foley et al. 
appear to have a normal Type la 


2010). These objects 
ight curve shape, ex¬ 


cept for the absence of a secondary maximum or “shoul¬ 
der” in near-IR bands. The reddest filter available for the 
HFFMTom light curve is F160W, which has an effective 
rest-frame wavelength of 6623A at z = 1.31, making it 
close to the rest-frame r band. Although the near-IR 
shoulder is more prominent in the rest-frame i band, we 
can see in Figured that the F160W light curve may sug¬ 
gest a weak or delayed near-IR shoulder. The second-to- 
last observation in F160W is ~ 1 a lower than predicted 
by the best-fit SALT2 and MLCS2k2 models. This is ten¬ 
uous evidence, but it would be consistent with a 2006bt- 
like light curve. 

In this case, however, the color of SN HFFMTom can 
rule out a classification as a 2006bt-like object. All 
members of this peculiar sub-class exhibit very red col¬ 
ors across all optical and near-IR bands, consistent with 
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a very cool photosphere (|Foley et al][2010j). Figure [7| 
shows that HFF14Tom is bluer than the SIN 2006bt pro¬ 
totype by at least 0.25 magnitudes in every color. The 
observed HFFMTom colors are fully consistent with a 
normal, blue Type la SN (the SALT2 model is shown). 


7.2. Possible Errors in Lens Modeling 

Having found no evidence for misinterpetations of the 
SN data, we are left to seek an explanation for the mild 
tension by scrutinizing the lens models. Here we evaluate 
four possible ways in which the input data or assump¬ 
tions of these lens models could lead to a systematic bias 
in the magnification. 


7.2.1. Unconstrained Mass Model Extrapolation 

As noted above, the position of SN HFFMTom just be¬ 
yond the edge of the strong-lensing region means that all 
models relying only on strong-lensing constraints must 
extrapolate the mass profile in order to make a predic¬ 
tion for the HFFMTom magnification. Let us suppose 
that the true cluster density profile for Abell 2744 hap¬ 
pens to have a sharp drop right at the edge of this core 
strong-lensing region. In that case these strong-lensing 
only models might well overestimate the mass interior to 
the HFFMTom position, and thus systematically overes¬ 
timate the magnification. 

One might expect that models incorporating weak- 
lensing constraints would be less susceptible to such bias 
from a sharply steepening mass profile. However, we 
find that collectively the four strong+weak lensing mod¬ 
els (shown as diamonds in Figurelb]) exhibit the same 
propensity to overestimate the magnification along this 
line of sight. This comparison therefore provides no evi¬ 
dence to support the extrapolation error hypothesis. 

With two versions of the CATS(vl) model, we also 
have a more direct test of the effect of introducing weak 
lensing constraints. The initial CATS(vl) model uses 
only strong-lensing constraints, and gives fi = 2.29+° ^, 
slightly higher than the measured value. A second itera¬ 
tion of this model, labeled here as CATS(vl.l), used the 
same strong-lensing features, but added in weak lensing 
constraints. The revised magnification of fj, = 2.62±0.18 
is further from the measured value, which serves to rein¬ 
force the conclusion that adding weak lensing constraints 
cannot, by itself, resolve this bias. 

Given the intermediate location of this SN, caught be¬ 
tween the strong-lensing and weak-lensing regimes, an 


flexion signal 

iGolclberg 

& Natarajan 2002 Goldberg 

& Bacon 20051 Irwin & 

Shmakoval 2006; Goldberg & 

Leonard 2007^ 

Flexion 

is a second-order lensing effect 


not lensed strongly enough to be warped into arcs or fil¬ 
aments. A flexion measurement near the SN HFFMTom 
sight line would be more sensitive to local gradients and 
substructure in the cluster mass profile. However, no 
such flexion constraints are yet available. 


7.2.2. Line of Sight Structure 

In addition to the dominant gravitational lensing 
from the foreground cluster, a distant source like SN 
HFFMTom will also be subject to the cumulative weak 
lensing effects due to uncorrelated large scale structure 


Weak Lensing” (Wong et al. |2011 Host 2012 

Collett 

et al. 2013; Greene et al. 2013 

Bayliss et al. 201 

1.41 |Mc- 

Gully et al. [2014 LbAloisio et a 

L 2014). This LSs effect 


which have a much longer path length over which to en¬ 
counter large scale structure and for highly magnified 
sources very near to the cluster lensing critical curve (be¬ 
cause LSS can perturb the position of the critical curve). 
For SN HFFMTom at the modest redshift of z = 1.35 
and far from the critical curve, LSS should increase the 
scatter in the len sing magnification at a level much less 
than 10% (D’Aloisio et al.[2014). The LSS effect is there¬ 


fore unable to account tor the observed discrepancy of 
~ 23% by itself. 

In addition, some sources of uncertainty are taken into 
account only by a subset of models. For example, the 
Wiliams model (using GRALE; Liesenborgs et al. 2006, 
2007; Mohammed et al. 2014) accounts for mass sheet 
degeneracy by introducing a free parameter in the recon¬ 
struction that represents an arbitrarily scaled mass sheet. 
The flexibility of the Williams model and the inclusion 
of this source of uncertainty contribute to its relatively 
large uncertainties compared to the other models. 


7.2.3. Nearby Cluster Member 

A component of the Abell 2744 lens that is particularly 
relevant to SN HFFMTom is the cluster member galaxy 
that lies just 578 north of the SN position (see Figure [I]). 
If the mass-to-light ratio or the profile of the dark mat¬ 
ter halo for this galaxy were significantly different from 
other cluster member galaxies, then its proximity to the 
SN sight-line might drive a bias in the magnification. We 
tested this hypothesis using the CATS(v2) model by al¬ 
lowing the mass of the nearby cluster member galaxy to 
vary as a free parameter in the model. We found that 
the change in the SN HFFMTom magnification predic¬ 
tion was less than A/i = 0.1. This additional dispersion 
is already included in the uncertainties quoted for that 
model in Table [5] and Figure [6] Furthermore, an erro¬ 
neous M/L value for the nearby cluster member galaxy 
could not explain the systematic shift of all lens models, 
because some do not incorporate cluster member galax¬ 
ies into their constraints at all (the free-form Bradac, 
Williams and Merten models). 


7.2.4. Cosmological Parameter Uncertainty 

Most lens models make a fixed assumption for the val¬ 
ues of cosmological parameters in a standard ACDM 
cosmology (typically: H m =0.3,flA=0.7,7Lo = f0 km s -1 
Mpc^ 1 ). This can introduce a systematic magnification 
error that is comparab le in magnitude to t hat of the sta¬ 
tistic al uncertainties (Zitrin et al. 2014 Bayliss et al. 
20151. Incorporating this cosmological parameter uncer¬ 


tainty would increase the model magnification errors, and 
therefore reduce the tension between models and obser¬ 
vati o n, but would no t resolve the overall systematic shift. 
Bayliss et al. (20151 find that for Abell 2744 the addi- 
tional cosmological uncertainty for magnifications /it ~ 2 
is < 10% which would still leave 6 models discrepant by 
at least 1.5cr. 


7.2.5. Misidentification of a multiple image 
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Jauzac et al. (2015) proposed a correction for the loca¬ 
tion ot i mage .'1.3. the third component of a tripl y-imaged 
galaxy ( Mer ten et al.||2011| at z — 3.98 ± 0.02 ^Johnson 
et al. |2014|). This possibly specious multiple image is 
only ~10'' from the SN HFF14Tom position, meaning a 
change in its location could have a large impact on the 
predicted magnification for SN HFF14Tom (see Johnson 
et al. 2014 for a quantitativ e discussion of this effect). 
However, Tauzac et al. (2015) also evaluated a version of 
the CATS(v2) model in which this single source is left in 
the original position and found that the location of image 
3.3 does not significantly affect the predicted magnifica¬ 
tions. Furthermore, dividing our 17 models into those 
using the original location vs those that adopt the new 
position, we find that both groups include models that 
are within 1 a of the observed /i as well as significantly 
discrepant points. 

7.2.6. Source Plane Minimization 

A general tendency towards high model magnification 
can result when lens models are minimized in the source 
plane (rather than the image plane). For a given model, 
modifying the surface mass density to have a shallower 
profile (declining less quickly with radius) will lead to 
an increase in the predicted magnifications across the 
field. As the magnification applies to both flux and area, 
the shallower profile also leads to a smaller source plane 
area, which necessarily brings delensed images of the 
same source closer together in the source plane. Thus, 
models that are naively optimized in the source plane 
will have a tendency toward shallower mass profiles and 
larger magnifications, which force all sources toward_the 
same location in the source plane (see, e.g., Kochanek 


1991 2006 Jee et al.]|2007[ Ponente & Diego 201l|). 


There are several strategies for mitigating this well 


avoid this bias entirely (Broadhurst et al. 

2005 Zitrin 

et al.||2009). The Sharon(v2), Zitrin-LTM, 

Zltrm-NFWT 


imization. Alternatively, one may guard against this 
bias by using a modified source plane minimization that 
is const rained not to gene r ate overly small de-lensed 


sources (Sendra et al. 2014 Diego et al. 2014). In our 
test set, a version of this second strategy is employed by 
the GLAFIC, Williams, and Diego models. For exam¬ 
ple, the Williams model uses source plane optimization 
in the GRALE algorithm, but avoids bias by maximizing 
the fractional overlap of extended source plane images. 
Finally, the effect of this bias should be reduced when one 
uses many strong lensing constraints over a wide range of 
redshifts. The CATS(v2) and CATS(v2.1) models have 
the largest number of strong-lensing constraints, using 
50 and 55 multiply-imaged systems, respectively. 

Although most of the tested models utilize one of these 
strategies to prevent a source plane minimization bias, 
Table [5] and Figure [6] show that these models do not 
uniformly deliver more accurate magnification predic¬ 
tions. The most direct test of this is in the Sharon 
and CATS model series. Both of these model series use 
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allows a choice ot how the minimization is done. For 
the Sharon series, the vl model used the faster source 
plane minimization approach, but the v2 model was con¬ 


structed with the more accurate image plane minimiza¬ 
tion. The CATS vl, 2, 2.1, and 2.2 models all used im¬ 
age plane optimization, while the vl.l employed source 
plane optimization (and also added in weak lensing con¬ 
straints). In neither of these model series do we see a 
uniform improvement in accuracy when optimization is 
done in the image plane. This test strongly suggests that 
the observed tension between the lens models and the SN 
observation cannot be entirely attributed to a bias aris¬ 
ing from source plane minimization. 

7.3. Quantity and Quality of Strong Lensing 
Constraints 

None of the possible errors in lens model inputs de¬ 
scribed above can completely account for the mild ten¬ 
sion between the observed and predicted magnifications. 
However, we do have several models that agree with 
the measured SN magnification. In this section we take 
a simple pragmatic approach and seek to identify any 
shared characteristics of the models that are most ac¬ 
curate for this particular sight-line. In this way SN 
HFFMTom may serve as a guide toward optimizing fu¬ 
ture lens models for magnification predictions, particu¬ 
larly in regions where strong-lensing constratins may be 
scarce. 

The addition of deeper imaging and new spectroscopic 
constraints from the HFF program and affiliated efforts 
might be expected to improve the accuracy of magni¬ 
fication predictions. However, in Figure [6] we saw no 
evidence for a universal improvement in model accuracy 
when moving from the pre- to the post-HFF models. In 
Figure[8]we attempt to refine this hypothesis by isolating 
the quantity and quality of strong lensing constraints. 

The left panel plots the tension between each model 
prediction and the measured SN magnification against 
the number of multiply-imaged galaxies used as strong- 
lensing constraints. The Bradac, Lam and CATSv2 
model series are highlighted here with connecting lines, 
as three examples where the number of systems used has 
changed substantially from one version to the next, but 
the model construction has remained fundamentally un¬ 
changed. There is no clear correlation here, and in fact 
for both the Lam and CATSv2 model series, we see the 
tension actually increases for a larger number of strong- 
lensing constraints. 

In the right panel of Figure[8]we examine an alternative 
explanation by plotting the model tension against the 
fraction of strong-lensing constraints that have a spec¬ 
troscopic redshift. For the SN HFFMTom magnification 
prediction, we see that all three of the testable model 
series (Bradac,Lam, CATSv2) become more accurate as 
the spec-z fraction increases]^] The spec-z fraction here 
serves as a crude stand-in for measuring how stringently 
the strong-lensing constraints are selected in these mod¬ 
els. 

The SN HFFMTom comparison therefore supports the 
notion that the quality of strong-lensing constraints is 
just as important as the quantity - at least for computing 
magnifications of sources outside the strong-lensing re- 

47 Of course this apparent correlation does not imply any di¬ 
rect causation. Consider the limit of just a single multiply-imaged 
galaxy constraint: the spec-z fraction could be 1.0 but the model 
would be wildly unconstrained. 
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Figure 8. The model tension as a function of strong-lensing constraints. The y axis in both panels marks the number of standard deviations 
in discrepancy between the measured SN magnification and the prediction from each lens model: (/i mo dei — MSN )/(<?■ model + <T 1 n) 1/2 - The 
left panel plots this tension against the total number of multiply imaged systems used, and the right panel plots the fraction or multiply- 
imaged systems that have a spectroscopic redshift constraint. Marker shapes and colors are the same as in Figure [6] In both panels, three 
lens model sequences are highlighted with connecting lines: the Bradac, Lam, and CATS(v2) models, which are the three examples where 
the number of strong lensing constraints and the fraction with spectroscopic redshifts changed substantially from one model version to the 
next. 


gion. Once a sufficient number of robust multiply-imaged 
systems have been included (perhaps ~10) most models 
will be dominated by syst emat ic uncertainties such as 
those described in Section |7.2| To improve the accu¬ 
racy of magnification predictions with a large number 
of strong-lensing constraints, one should employ a strict 
vetting of the multiply-imaged systems. Each system 
used should ideally have a spectroscopic redshift, or at 
least a very well constrained photometric redshift. The 
primary value in this strict vetting is to avoid erroneous 
redshifts being assigned to multiply-imaged systems that 
are unconstrained in redshift space. Such redshift er¬ 
rors are especially problematic for relatively low -2 back¬ 
ground systems: at lower redshifts the model-derived de¬ 
flection angle and magnification show a more rapid vari¬ 
ation with the redshift of the background galaxy than 
for similarly positioned high redshift galaxies. This con¬ 
clusion should be unsurprising, as for example, |Bayliss| 


et al. (2015) have shown that cosmological uncertain¬ 


ties are smaller for models that use a large fraction of 
stro ng-lensing constraint s with spectroscopic redshifts, 
and Johnson et al. (20141 have shown a case study where 
adding"spcctroscopic redshift information changes the in¬ 
ferred magnification by ~10%. 

It is relevant to note that the three model series high¬ 
lighted here (Bradac, CATSv2 and Lam) also happen to 
be the three models that are not truly “blind,” in that 
the modelers were aware of preliminary analysis of the SN 
magnification when constructing the latest version. How¬ 
ever, none of these models used the SN magnification as 
an input constraint, and none were consciously tuned to 
match the SN magnification. In fact, for the Bradac(v2) 
and Lam(v2) cases, the modelers did not have access to 
a final measurement of the SN magnification at the time 


their models were generated. 

8. SUMMARY AND CONCLUSIONS 

The appearance of a Type la SN behind a massive 
galaxy cluster provides a rare opportunity to use a stan¬ 
dard candle for a direct measurement of the absolute 
magnification due to gravitational lensing. The discov¬ 
ery of SN HFF14Tom in the HFF imaging of Abell 2744 
offers the first chance to apply this test on a cluster with 
multiple publicly-available lens models. We have found 
that the spectrum and light curve of SN HFF14Tom 
are well matched by templates of a normal Type la SN 
at 2 = 1.3457. Using the two most prevalent Type 
la SN light curve fitters, SALT2 and MLCS2k2, we 
get a consistent measurement of the distance modulus 
(Table [3]) . Using a cosmology-independent comparison 
against a sample of unlensed Type la SNe at similar red¬ 
shifts, we find that SN HFF14Tom is ~ 0.7 magnitudes 
brighter than the field sample would predict. Attributing 
this difference to the gravitational lensing magnification 
(and accounting for the intrinsic scatter in luminosity of 
the Type la SN population), we have derived a consis¬ 
tent measured magnification of \x salt 2 = 1-99 ± 0.38. 
A*MLCS2k2 = 2.03 ± 0.29, from the two light curve fitters. 

Taking advantage of the availability of 17 well- 
constrained lens models for the Abell 2744 cluster, we 
have used SN HFF14Tom to ask how accurately these 
lens models can predict the magnification along this line 
of sight. We find that these models are consistent, and 
fairly accurate, collectively predicting /j = 2.5 ± 0.4, 
within ler of the measured value. This is encourag¬ 
ing, and reinforces the quality and value of these public 
lens models for studying magnified background objects. 
However, we note that all but two of the tested models 
provide a median magnification prediction that is larger 
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than the measurement from the SN, and all the mod¬ 
els with significant discrepancy are too high. This may 
be an indication that there is a small systematic bias at 
work, at least in some of the models. 

We have speculated on what could be the origin for 
such a systematic bias, first considering and rejecting 
three possible explanations that presume an error in the 
interpretation of the available SN data. 

1. Redshift: We reject the possibility that a redshift 
error is the primary cause of the discrepancy, as the 
redshift evidence is well supported by multiple lines 
of evidence from both the SN and the presumed 
host galaxy. Changing the redshift within the con¬ 
straints of these complementary observations does 
not resolve the tension. 

2. Dust: We find it implausible that there is sufficient 
dust in the cluster or elsewhere along the line of 
sight to account for the magnification discrepancy. 
A dust-induced bias is also disfavored by the very 
blue color of SN HFF14Tom. 

3. Misclassification: The combination of spectro¬ 
scopic and photometric evidence strongly supports 
our classification of SN HFF14Tom as a normal 
Type la SN. The most plausible mis-classffication 
would be that the object is a peculiar Type la of the 
SN 2006bt-like sub-class. Although the light curve 
shape could allow this possibility, the blue color of 
HFF14Tom can once again reject this alternative. 

Turning to the lens models, we have considered 6 ways 
in which a bias might be introduced into the lens models: 

1. Mass profile extrapolation errors 

2. Large scale structure along the line of sight 

3. Peculiar mass for a nearby cluster member 

4. Cosmological parameter uncertainty 

5. Misidentification of a multiple image 

6 . Source plane minimization 

Although each of these could plausibly introduce a small 
systematic bias, we find that none are likely to be signif¬ 
icant enough (or universal enough) to completely resolve 
the tension between models and observations. 

Finally, we have examined whether there is a simple 
prescription for the kind of strong lensing constraints 
that are most likely to deliver an accurate magnification 
prediction for this sight-line. We find that the number of 
multiply imaged systems used is not in itself predictive 
of model accuracy. Rather, it is the fraction of multiply 
imaged systems that have spectroscopic redshifts that is 
most correlated with model accuracy. A reasonable in¬ 
terpretation - at least for this cluster and this particular 
set of models - is that the spectroscopic redshift frac¬ 
tion serves as an effective proxy for the “quality” of the 
strong lensing constraints. This quality of the input data 
appears to be a key ingredient for deriving accurate mag¬ 
nifications outside the strong-lensing regime. 


We have evaluated here only a single object behind 
a single cluster, with a minor tension between the ob¬ 
served magnification and the model predictions. This is 
not in and of itself a cause for alarm. Previous analy¬ 
ses of lensed Type la SN found no significant discrep¬ 
ancy between the observed Type la SN magnifications 
and t he predictions from lens models (P14 Nordin et al. 
20141, albeit with a much smaller set of lens models be- 


ing tested. The observed systematic bias for HFF14Tom 
is small, and many of the lens models that deviate 
from the measured magnification are preliminary models 
that have not been updated to include all of the HFF 
data. Future revisions of the lens models for Abell 2744 
could either incorporate the observed magnification of 
HFF14Tom as a new model constraint, or can revisit 
this test to evaluate whether the bias persists. 

The general problem of combining predictions from 


ing et al. 1999 Liddle 20091. In a recen 


Press 1997 
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; example, 
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et al. 12013) examined the photometric redshifts from 


multiple galaxy SED fitting codes using the same input 
data. They found that the photo-z accuracy and pre¬ 
cision can be substantially improved by combining the 
outputs using a sum of probability distributions, a hier¬ 
archical Bayesian method, or even a technique as simple 
as taking the straight median from all available models. 
It may be that lensing magnifications are similar, in that 
they can in general be improved through a simple combi¬ 
nation of lens model predictions. Such model averaging 
will not, however, resolve a systematic bias that is shared 
by many models. 

8 .1. Future Work 

A promising avenue for exploring the origins of such 
systematic biases in cluster lens models is through the 
use of simulated lensing data. One can start with very 
deep high-resolution multi-band imaging on an unlensed 
field that has a fairly complete spectroscopic redshift 
catalog, such as the Hubble Ultra Deep Field. Then a 
simulated galaxy cluster is placed in the field, and the 
background galaxies are distorted into arcs and multi¬ 
ple images using a well-defined lensing prescription. The 
artificially lensed images can then be distributed to lens 
modeling teams who attempt to reconstruct the (known) 
mass profile of the simulated cluster. This exercise has 
recently been pursued with a set of synthetic clusters sim¬ 
ilar to those observed in the HFF program (Meneghetti 
et al. in prep), u sing techniques similar to Meneghetti] 
et al. (2010 2014). Preliminary analysis of this simula- 


tion comparison suggests that magnifications can be sys¬ 
tematically overestimated for sources that lie outside the 
strong-lensing region along the minor axis direction of a 
simulated cluster that is very elongated, similar to Abell 
2744. This indicates that some lens models significantly 
underestimate the ellipticity of the mass distribution on 
large scales for that simulated cluster, which is consistent 
with a recent lens model compa rison using 25 rea l clus¬ 
ters from the CLASH program (JZitrin et al.]|2015). Our 
analysis of SN HFF14Tom here indicates that such sim¬ 
ulation efforts and head-to-head lens model comparisons 
will be an important step for moving toward precision 
science with cluster-lensed sources. 

Increasing the sample of SNe behind clusters like those 
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in the HFF program with rich lensing constraints and 
deep imaging - would allow the test described here to 
be repeated and refined. With 10 or 100 such objects, 
it would be possible to see whether the SN HFF14Tom 
H discrepancy is simply an outlier, or an indication of 
a more pernicious systematic error. Any cluster that 
has been vetted by pencil-beam magnification tests using 
lensed Type la SN will be able to provide a more reliable 
measure of the magnifications for very high redshift ob¬ 
jects. Similarly, a broad sample of SNe like HFF14Tom 
would help to define the preferred lens modeling method¬ 
ology by highlighting any models that consistently per¬ 
form well in Type la SN lensing tests. The ongoing 
FrontierSN program will discover and follow any more 
highly magnified SNe that appear behind the Frontier 
Field clusters. Unfortunately, the HFF survey is not 
designed with high-z transient discovery as a primary 
science goal, so the FrontierSN effort will likely add no 
more than 1-3 new lensed Type la SN. Further imaging 
of strong-lensing clusters with HST or the James Webb 
Space Telescope (JWST) could enable a larger sample 
to be collected, especially if the filters and cadence are 
optimized for detection of Type la SN at z > 1. Massive 
clusters such as Abell 2744 will continue to be attrac¬ 
tive as cosmic telescopes, allowing the next generation of 
telescopes to reach the faintest objects in the very early 
universe. The puzzling bias revealed by SN HFF14Tom 
supports a concerted effort to improve these lenses with 
further examination of lensing systematics through sim¬ 
ulations and collection of a larger sample of magnified 
SNe. 
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